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Given a polynomial system /, a fundamental question is to determine if / has real roots, 
ly-^ Many algorithms involving the use of infinitesimal deformations have been proposed to 

answer this question. In this article, we transform an approach of Rouillier, Roy, and Safey 

El Din, which is based on a classical optimization approach of Seidenberg, to develop a 
I I homotopy based approach for computing at least one point on each connected component 

of a real algebraic set. Examples are presented demonstrating the effectiveness of this 

parallelizable homotopy based approach. 
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CN 1 Introduction 

> 

Computing real roots of a polynomial system is a difBcult and extremely important problem. 
In many applications in science, engineering, and economics, the real roots are the only ones of 
interest. Due to the importance of this problem, many approaches have been proposed. Two 
approaches are the cylindrical algebraic decomposition algorithm [T7 and so-called critical point 
methods, such as Seidenberg's approach of computing critical points of the distance function 
[41j . The cylindrical algebraic decomposition algorithm has doubly exponential complexity in 
the number of variables. However, using the idea of Seidenberg and related ideas developed in 
[5l[T6l[2Tll22ll2Il[39], algorithms with asymptotically optimal complexity estimates for computing 
at least one real point on each connected component of a real algebraic set were developed. Other 
related approaches for computing real roots are presented in [TJ [31 131 SO] and the references 
therein. The approach presented here will transform the algorithms presented in [TJ 20] into a 
homotopy based algorithm. 

Several homotopy based algorithms have been proposed to compute real roots of a polynomial 
system. The algorithms in [32 and 13_ utilize critical point methods to decompose the real 
points of a complex curve and a complex surface with finitely many singularities, respectively. An 
algorithm for directly computing only the real roots that are isolated over the complex numbers 
is presented in 12J . The complexity of this approach depends upon the fewnomial structure of 
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the given polynomial system. The approach presented below is not restricted to low-dimensional 
cases and the real roots are not assumed to be isolated over the complex numbers. 

Two other nonhomotopy based algorithms are presented in [30 and [IS]- The approach 
in [30] (see also [3T]) uses semidefinite programming for computing real roots. This algorithm 
computes every real root assuming the number of real roots is finite. The approach in [151 uses 
tools related to maximum likelihood estimation in statistics for computing real positive roots of 
certain types of polynomial systems. 

The rest of the article is structured as follows. The remainder of this section describes the 
needed concepts from complex, real, and numerical algebraic geometry and a brief introduction 
to Puiseux series. Section [2] describes the homotopy based approach with examples demonstrat- 
ing the algorithm in Section |3] 

1.1 Algebraic sets and genericity 

Let / : C" be a polynomial system and V(/) = {x e | f{x) = 0}. The set V(/) C 

is called the algebraic set associated to f. A set X C is called an algebraic set if there exists 
a polynomial system g : — C" such that X — V{g). An algebraic set X C is reducible if 
there exists algebraic sets Y,Zc , which are proper subsets of X, such that X ^ YU Z. An 
algebraic set is irreducible if it is not reducible. For an irreducible algebraic set X, the subset 
of manifold points X^-^g is dense in X, open, and connected. The dimension of an irreducible 
algebraic set X is the dimension of X^^g as a complex manifold. 

On irreducible algebraic sets, we can define the notion of genericity. 

Definition 1 Let X C be an irreducible algebraic set. A property P is said to hold 
generically on X if the subset of points in X which do not satisfy P are contained in a proper 
algebraic subset of X. That is, there is a nonempty Zariski open subset U oi X such that P 
holds at every point in U. Each point in U is called a generic point of X with respect to P. 

Since every proper algebraic subset of C is a finite set, a property P holds generically on C 
if P holds at all but finitely many points in C. 

Every algebraic set X can be written uniquely (up to reordering) as the finite union of 
inclusion maximal irreducible algebraic sets, called the irreducible decomposition of X. That is, 
there are irreducible algebraic sets Ai, . . . ,Ak such that 

fc 

X ^[J Ai and Ai (/_ Aj for i ^ j- 

i=l 

Each Ai is called an irreducible component of X. 

The dimension of an algebraic set is the maximum dimension of its irreducible components. 
An algebraic set is called pure- dimensional if each irreducible component has the same dimension. 
The pure i-dimensional component of an algebraic set is the union of the irreducible components 
of dimension i. In summary, the algebraic set V(/) has an irreducible decomposition of the form 

dimV(/) dimV(/) fc, 

v(/)= U U U^^^^ (1) 

1=0 i=0 j = l 

where Vi is the pure z-dimensional component of V(/) and each Vij is a distinct i-dimensional 
irreducible component. 
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1.2 Decomposition of real algebraic sets 

A real algebraic set are subsets of which arise as the intersection of algebraic sets in 
with M^. That is, a set X C is a real algebraic set if there is an algebraic set Y C such 
that X = Y C\ M^. For a polynomial system / : — >■ K", the real algebraic set associated to f 
is Vk(/) = V(/) n = {x e I fix) = 0}. 

Consider the algebraic set X = V{y^ — {x — \)) <Z . It is easy to verify that X is an 
irreducible algebraic set and, hence, both X and X^og are connected. However, the real algebraic 
set n is not connected. This example suggests that we should consider decomposing real 
algebraic sets into connected components. 

A real algebraic set X C can be written uniquely (up to reordering) as the disjoint 
union of finitely many path-connected sets Ci , . . . , C such that Ci and V \Ci are both 
closed in the Euclidean topology on K^. Each Ci is called a connected component of X and 
one can verify that it is a semi- algebraic set. Expanded details regarding real algebraic sets and 
decompositions can be found in [U [H] . 

To demonstrate this decomposition and contrast it with the irreducible decomposition of 
algebraic sets, consider the algebraic sets X — V{y^ — a;^(x — 1)), Y = V{x — y), and Z — XUY 
with corresponding real algebraic sets — X O M.'^ , Yg^ — Y O M."^ , and Zr ^ Z OM.'^. It is easy 
to verify that X and Y are irreducible algebraic sets with Z clearly being a reducible algebraic 
set. The set consists of two connected components, namely Ci = {(0,0)} and a connected 
curve C2 = yR\Ci. Since the real algebraic sets Yr and Zr are connected, 1r and Zr each have 
only one connected component. 



1.3 Puiseux series 



Since we will utilize Puiseux series in Section [2] we will provide a brief review here. For more 
detailed information, see |3]. 

The field of algebraic Puiseux series over C is 



. j>jo 



J/l 



jo G Z, g e N, e C with Oj^ 7^ 



To simplify the notation, we shall define aj — for all j < jo- An element in C(e) is bounded if 
jo > and infinitesimal if jo > 0. The subset consisting of bounded elements, denoted Cb{e), is 
a ring which is naturally mapped to C by the ring homomorphism limo defined by 

limo ^ o.je^^'^ = oq. 

j>jo 



1.4 Numerical irreducible decomposition and witness sets 

Let / : C" be a polynomial system. A numerical irreducible decomposition of V(/), 

first presented in [48 , is a numerical decomposition analogous to ([T]) using witness sets (see [49l 
Chaps. 12-15] for more expanded details). Suppose that V is the pure i-dimensional component 
of V(/) with d = degV". For a fixed generic i-codimensional linear space H C C^, wc have 
that V O H consists of d points. Let L : — C* be a system of linear polynomials such 
that V{L) = H. The set V n V{L) = V fMI is called a witness point set for V with the triple 
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yy = {/, L,V D V{L)} called a witness set for V. A numerical irreducible decomposition of V(/) 
is of the form 

dimV(/) ki 

U Uw^^ (2) 

i=0 3 = 1 

where Wi.j is a witness set for a distinct i-dimensional irreducible component of V(/). We 
note that the union of witness sets in ([2| should be considered as a formal union. Numerical 
irreducible decompositions can be computed using the algorithms presented in [6l [26l |42j |44] |45| 

ssmziiisj. 



1.5 Trackable paths 

Numerical homotopy methods rely on the ability to construct homotopies with solution paths 
that are trackable. The following is the definition of a trackable solution path starting at a 
nonsingular point adapted from [25 . 

Definition 2 Let H{x, t) : x C he polynomial in X and complex analytic in t and let 
X* be a nonsingular isolated solution of H{x, 1) = 0. We say that x* is trackable for t £ (0, 1] 
from t = 1 tot = using H{x, t) if there is a smooth map ^x* '■ (0, 1] such that ^j;. (1) = x* 

and, for t G (0, 1], £,x*{t) is a nonsingular isolated solution of H{x,t) — 0. 

The solution path starting at x* is said to converge if limf_^o+ (t) G C^, where limt_^Q+ S^x* {t) 
is called the endpoint (or limit point) of the path. 



2 Real points on an algebraic set 

Let / : — > M" be a polynomial system and V C Vc(/) be a pure d-dimensional algebraic 
set. The main problem we consider is, given a witness set {f,C, W} for V, compute a finite set 
of points which contains at least one point on each connected component of V]r(/) contained in 
V. We note that if c? = 0, then V — W so that one can compute the real points in V simply by 
considering the finitely many points in W. Hence, we will assume that d > 0. 

We will also reduce to the case n = N — d. One way to always reduce down to this case is 
to consider the polynomial g which is the sum of squares of /, that is, 3 = /i + ■ ■ • + with 
Vg = V{g). If T is a finite set of points which contains at least one point on each connected 
component of V]R(g), then T nV contains at least one point on each connected component of 
"^'r(/) = Vr{9) contained in V. The set Tn ^ can be computed from T and a witness set for V 
using the homotopy membership test j45| . 

We summarize the assumptions in the following statement. 

Assumption 3 Let N > d > 0, f : M^^'' be a polynomial system, and V C V(/) be a 

pure d-dimensional algebraic set with witness set {/, £, W}. 

The following lemma considers the solutions of f{x) = z for z e C^^'*. 

Lemma 4 With Assumption^^ there is a nonempty Zariski open set Z C C^^'* such that, for 
every z (z Z , V(/ — z) is a smooth algebraic set of dimension d. 
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Proof. Let r = diin/(C^) and c = N — r, which are caUed rank of f and the corank of f 
respectively [Ml §13.4]. Since /(C^) C C^"'*, we have r < N - d and hence d < N - r = c. 
Since V(/) has a component of dimension d, Theorem 13.4.2 of H3] yields that d> c. Therefore, 
c — d and r = N — d. The lemma now follows immediately from Lemma 13.4.1 of |49j . □ 

Lemma |4] permits the use of continuation techniques as stated in the following theorem. 

Theorem 5 Suppose that Assumption^ holds. Let z G M^^'',, 7 G C, y G \ V]r(/), a G 
C^-''+S and H -.C^ X C^-'^+i x C ^ c2W-d+i ^^^^ homotopy defined by 



H{x,X,t) 



f{x) - t-/z 

Ao(a; -y) + AiV/i(a:)^ + • ■ ■ + Xn-ci"^ fN-d{x f 
aoXo + h a^^dXN-d - 1 



(3) 



where f{x) — [fi{x), . . . , fN-d{x)]'^ such that following statements hold. 

1. The set S (Z xC^^"^^^ of roots of H{x, X,l) is finite and each is a nonsingular solution 
ofH{x,X,l)=0. 

2. The number of points in S is equal to the maximum number of isolated solutions of 
H{x,X, 1) — as z, 7, y, and a vary over sets C^^"^, C, , and C^^'^'^^, respectively. 

3. The solution paths defined by H starting, with t — 1, at the points in S are trackable. 

4. IfTT{x,X) = X, 

E = {lim(_^o+ '?s(i) I s Cz S and the solution path converges} , and 
El = {lim(_^o+ ""(^slO) I s G 5* and the path 7r(^s) converges} , 

we have Ei = Tr{E). 

Then, Ei OV D M.^ contains a point on each connected component 0/ V]r(/) contained in V. 

The homotopy H defined in ([3| is based on the classical approach of Seidenberg [IT]. If 
y G \ Vr(/), consider the quadratic polynomial 

N 

dy(x) = (x - y)^(x -y) = ^(x, - y^)^ 

i=l 

and the optimization problem 

(P) mm {dy{x) I XG Vk(/)}. 

We want to compute the points on V(/) for which Wdy{x) — 2{x — y)^ and Wg{x) are linearly 
dependent. The approach in [40j for hypersurfaces uses determinants to describe this linear 
dependence condition, while the approach in Theorem [5] uses auxiliary variables A. In particular, 
the polynomial system Gf^y : x P^"'' C^^"'' defined by 



Gf^y{x,X) 



Xoix -y) + XiVfi{x)^ + ■■■ + XN-d'^fN-dixy 



(4) 



comprises the Fritz John conditions 29; for problem (P) and provides necessary conditions for 
optimality. That is, if ^ is a local minimizer for problem (P), then there exists A G P^^*^ such 
that (C,A) G ViGf^y). 
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Clearly, x E V(/) such that 

mnk[x-y V h(x)'^ ••• V fN-d{xf ] < N ~ d 

if and only if there exists A G P^~'' such that (x, A) G '^{Qf,y)- A point x £ T^(V{Gf.y)) is called 
a critical point of the distance function with respect to /, where 7r(a;, A) = x. 

Consider Sing(/) = {x G | rank Jf{x) < N — d}, where J f{x) is the Jacobian matrix of 
/ evaluated at x. If Sing(/) is positive dimensional, then V{Qf,y) is also positive dimensional. 
By using Lemma [4j we can consider smooth algebraic sets thereby allowing the computation of 
finitely many points in Vk(/) containing the points of interest. 

The following lemma will be used to complete the proof of Theorem [5] 

Lemma 6 Suppose that Assumption^^ holds. Let e be an infinitesimal, y G M.^ \ Vr(/), z G 
jjiV-d ^j^/j =^ 0, and fe{x) = f{x) — ez be such that V{Qf^^y) is finite and \V{Gf^.y)\ is equal 
to the maximum number of isolated solutions as y and z varies over the sets and C^^'', 
respectively. Then, 

1. -l/climo(V(/e)nCf,(e)^), 

2. limo (V(/,) n Cb(e)^) n = limo {V{fl - e^zj, • . . , /^.^ - e^z%_^) D Rb{e)^), and 

3. limo nC{,(e)^) n contains a point in each connected component ofV-^{f) 
contained in V where 7r(a;, A) = x. 

Proof. This setup implies that V(/e) is a c?-dimensional smooth algebraic set for which we 
clearly have V C hmo(V(/e) n Cb(e)^) yielding Iteni[l] Item [2] follows from the fact that G M 
and 

limo (V(.9 - e) n Cb(e)^) n = limo (V(5' - e') H Mb(e)^) 

for any polynomial g : M. Item [s] follows by using the same proof as Lemma 3.7 in [JO] 

with the replacement of Lemma 3.6 of [10] with Items [T] and [2j □ 

Before we prove Theorem [sj we note that the polynomial system Qg,y{x,\) defined in Q, 
has A G P^"''. The polynomial system ff(a:,A,0) defined in ([3| has A G C^~''+^ restricted to 
the Euclidean patch defined by OqAo + • • ■ + C(N~d^N-d — 1- Item [3] in Theorem [5] enforces 
that this Euclidean patch is in general position with respect to the finitely many solution paths. 
Therefore, we can use the results of Lemma |6] in the following proof of Theorem [5] 

Proof of Theorem [5} Let e be an infinitesimal and fe — f — <^z. Item [2] yields that \S\ = 
\V{Qf,,y)\ < 00. The result will follow from Lemma[6]upon showing 

£^1 =limo(7r(V(^/^,,))nCfc(e)^). (5) 

We will deduce ([s]) by comparing the polynomial systems Gf,,y and 



f{x) — ez 

Ao(a; ~y) + XiV fi{xY + ■■■ + XN-d^fN-dix^ 
aoAo H 1- aN-d^N-d - 1 



Item[2 also yields that \S\ = \V{Gf^^y)\ = \V{Gf^ y)\- In particular, by abuse of notation regarding 
TT, wenave 

AViGf.^y)) - HnSty))- (6) 



6 



Since there arc finitely many honiotopy paths, there exists < to < 1 such that all of the 
honiotopy paths for H with < t < 2to are described by the points in V{Gf y) C C(e)^^~'' 
by replacing e with This yields that the set of limit points of the homotopy HQ(x^X,t) :— 
H{x, A, {1 — t) ■ to) starting at the roots of H{x, A, to) is 

r = limo (V(C/;^,,)nCb(e)^+^). 

Since, by Items [T] and [sj the homotopy paths of H are nonsingular for t £ (0,1], coefficient- 
parameter continuation [35] yields that T = E. Item |4] yields 

^{E) = IT (hmo {V{gi,y) n Cfc(e)^+2)) = limo {A^iGL^y)) n Cb(e)^) = E^. 

This equation together with ^ yields ([5|. □ 

We note that in the hypersurface case, that is n = N — d — 1, if / has degree 2k, the 
2-homogeneous Bezout count yields that 

IS"! < K{N, 2k) ■=N ■2k- {2k - 1)^-^ 

In particular, Vr(/) can have at most K{N,2k) connected components and hence K{N,2k) 
bounds the number of real roots of / that are isolated over M. This bound is only N times 
larger than the bound obtained in [TH Prop. 11.5.2]. 



2.1 An algorithm 

Theorem [5] yields an approach for computing a point on each connected component of Vr(/). Be- 
fore presenting an algorithm which implements the ideas of this theorem, we state two remarks. 
First, Item[2]of Theorem [5] holds for a nonempty Zariski open set of C^"'^ x C x x C^"'*+^ 
The following algorithm assumes that the given point {z, 7, y, a) lies in this Zariski open set. As 
part of the procedure, it computationally verifies Items [l| [3] and [4] of Theorem [5] hold. Second, 
the use of 7 is based on the "Gamma Trick" [491 Lemma 7.1.3] first introduced by Morgan and 
Sommese |34j . 

Second, since there exist many suitable methods to compute the start points S, the following 
algorithm does not directly specify which one to utilize. Nonetheless, to improve efficiency in 
this computation, the method should, in some way, utilize the natural 2-homogeneous structure. 

Procedure [v,R] = RealPoints(/, W, z, 7, y, a) 

Input A polynomial system / : M^"'', a witness set W for a pure d-dimensional algebraic 

set V C V(/), z e M^-'', 7 e C, y e \ Vr(/), and a G C^-'^+i such that Item [2] of 
Theorem [5] holds. 

Output A boolean v which is true if Items [l] [3j and[4]in Theorem [5] have been computationally 
verified, otherwise false. If v is true, i? is a finite subset of M containing a point on each 
connected component of the real algebraic set Vr(/) contained in V. 

Begin 

1. Construct the homotopy H defined in (|3|. 

2. Compute the solutions S of H{x, A, 1) — 0. 
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(a) Use S to verify that Item [T] of Theorem [5] holds. If it does not hold, Re- 
turn [false, ^. 

3. Track the solution paths of H starting at each point in S to compute the sets E and 
El defined in Theorem [5) 

(a) If the tracking fails for a path or Tr{E) ^ Ei where 7r(x, A) = x, Return [false, 0] . 

4. Use the homotopy membership test to compute the set R consisting of the points in 
El n contained in V. 

Return [true,R]. 

Since the endpoints E computed in Step [S] may be singular solutions of H{x,X,0) = 0, the 
use of an endgame, e.g., [51 [55J [371 , together with adaptive precision tracking [71[S1[TT] may 
be required to accurately compute them. Also, Step [s] should use the method of [33] to avoid 
infinite length paths. 

Example 7 To illustrate the algorithm for a hypersurface, consider the polynomial f{xi,X2,X3) ~ 
{xi + xs)^ + x\ with V = V{f). Clearly, Vm(/) = {{a, 0, -a) | a G M} C Sing(/). Item[2]holds 
with z = 1, 7 = 2 + 3i, 



3/8 
5/9 
1/3 



and a 



1/2- 
6/7 H 



2i/3 



where i 



Let H be the homotopy defined by (|3|. 

• ForStep[2) we used a standard 2-homogeneous homotopy, which required tracking K{3, 2) = 
6 paths, to compute the set S consisting of the four nonsingular solutions of H{x, A, 1) — 0. 

• The four paths tracked in Step[3j which started at the points in S, all converged with the 
endpoints of the two paths ending at the real point coinciding. In particular, E and Ei 
both consist of three points with Tr{E) = Ei where 



E^ = 



(1/48, 0,-1/48), (-1/3 + 5i/9, 10/9 + 17i/2A, -3/8 + 5i/9), 
(-1/3 - 5i/9, 10/9 - 17i/24, -3/8 - 5i/9) 



^ = {(1/48,0,-1/48)}. 

-1/48) is the minimizer of the distance between the 



• Since V = V(/), we have R = Ein] 

It is easy to verify that the point (1/48, 
point y and V]^{f), as shown in Figure [ij 

Example 8 To illustrate the algorithm for an algebraic set, consider the polynomial system 



/(^) 



92{x) 



■r29z{x) 



where g{xi,X2,xz) = 



{xl - X2){xl 

{xiX2 

{xixs, 



X3)(a;2 

X2){Xi 



l){xi^ 

-1)(^2 



2) 
3) 



with ri = 1/3 and r2 = 1/7. We want to investigate real points of the cubic curve V = 
{{xi,Xi,xl) I a;i G C} C V(/) which is input into RealPoints via a witness set W. Item [2] 
holds with 

\fg , 7 = 3/ll-z/13, 
Let H be the homotopy defined by ([s]). 



1/4 ■ 




1/3- 


-z/7 






1/6 


, and a = 


6/lH 


-3i/4 


where i = \ 




-3/2 




2/3- 


7i/8 
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Figure 1: Plot of Vr(/) and the point minimizing the distance between y and V]i{(/) 

• We used a standard 2-homogeneous homotopy, which required tracking 300 paths, to 
compute the set S in Step [2] consisting of the 95 nonsingular solutions of H{x, A, 1) = 0. 

• All 95 of the paths tracked in Step [3] starting from the points in S converged with the set 
El n consisting of 15 points. 

• The homotopy membership test yields that 7 of these 15 points lie on V . 

The point of minimum distance on y n to y is approximately (0.168, 0.028, 0.005), which is 
displayed in Figure [2] with the other 6 points on V . 




Figure 2: Plot of F n M and the point minimizing the distance between y and y n M 
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3 Examples 



The following examples were run using the software package Bertini vl.3.1 |10j on a server 
having four 2.3 GHz Opteron 6176 processors and 64 GB of memory that runs 64-bit Linux. 
The serial examples used one core while the parallel examples used one manager and 47 working 
cores. For the nonsingular solutions, we utilized alphaCertified [231 [M] to certify reality. For 
the singular solutions, we determined reality based upon the size of the imaginary parts using 
two different numerical approximations of the point. 

3.1 Hypersurface example 

Consider the polynomial provided in Example 5 of |40| . namely 

fiu2,U3,U4, U5) = IIOU5W4U3 + I9OU5U4U3 + 8OU4U3 + 8OU5U3 + 27OU5U4U3 + 160u|m3 
-I-8OM5U3 + 80m4u| - 32m4U§M2 - 32it3M2 - 80u|w2 - I28M5U4U2 

— I6OW5M3U2 — II2W4U3U2 — 64u|u| — 8OU5U2 — 32U3W2 + 60u|m4 
+220u5ul + 160u| + 67M5U4W3 + 136u|u3 - 24u5U§ - 88u4Mi - 64u| 

— IOOW5M2 + 32M5M4U2 + 96M4U2 — 228W5U3U2 — IO8M4U3M2 — 120u|u2 
-I-2OM5U2 + 96W4W2 — 56M3U2 + IIOM5U4 + 80ul + 48U4U3 — 32^3 
-|-30M5tt2 + 48U4U2 — 2OM3U2. 

The approach of [ID] computes 26 real points on the hypersurface which contains at least one 
point on each connected component of Vr(/) using y — 0. Since y ~ does not satisfy the 
hypotheses of Theorem^ we used y = [ 4/3 —9/5 —5/7 8/9 ] to compute at least one 
point on each connected component of Vi[{(/). In particular, we used serial processing with 
RealPoints taking V = V(/), z = 1, and 7 G C and a G to be random of unit length. 
Let H be the homotopy defined by 

• For Step [2] we used a 2-homogeneous regeneration [25] to compute the set S consisting of 
the 151 nonsingular solutions of H{x, A, 1) = in 10 seconds. 

• For Step [3] each of the 151 paths converged with the set Ei D consisting of 28 distinct 
points, which was computed in 4 seconds. 

• Since V V{f), R ^ Ei nR^ which consists of 28 points. 

Wc note that since |Sing(/)| < 00, we could directly compute V{H{x, A, 0)) using a standard 
2-homogeneous homotopy, which requires the tracking of 432 paths. Bertini performed this 
computation in serial in 27 seconds which yielded the same set R of 28 real critical points, as 
required by theory [35) . 

3.2 An example from filter banks 

Consider the polynomial system named F633 [19| that was considered in [1], which is available 
at [20]. This polynomial system consists of 9 polynomials in 10 variables. Since two of the 
polynomials are linear and linearly independent, we utilized intrinsic coordinates to reduce the 
number of variables to 8 and the number of polynomials to 7, all of which are bilinear. Since 
these 7 polynomials are not independent, we further reduced down to a system of 6 bilinear 
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polynomials in 8 variables, namely 



where 

g{xi, . . . ,X4,2/i, . . . ,2/4) = &ixiy2+xiy3 + x2y3 + xiy4: + X2y4 + X3y4)+'iixi+x2 + xa + x4) + 13. 

The algebraic set V(/) is an irreducible surface of degree 32. We used RealPoints to compute 
a set of points containing a point from each connected component of Ve(/) by taking 

y=[l/5 -3/4 -2/3 7/9 -4/7 12/13 1/2 -10/11 

z G M^, 7 € C, and a G to be random of unit length, and W a witness set for V — V(/). Let 
H be the homotopy defined by 

• For Step [2] we used a standard 2-homogeneous homotopy, which required tracking 1792 
paths, to compute the set S consisting of the 274 nonsingular solutions of H{x, A, 1) = 0. 
This computation took 120 seconds in serial (4 seconds in parallel). 

• For Step [s] each of the 274 paths converged with the set Ei n consisting of 36 distinct 
points. This computation took one second in serial. 

• Since V — V(/), R — EiD which consists of 36 points. 

In Step [2) we could have used a 3- homogeneous homotopy since the system itself is naturally 
2-homogeneous. However, this would increase the number of paths from 1792 to 1960. Also, 
since Sing(/) — 0, we could directly compute V{H{x, X,0)) using a standard 2-homogeneous 
homotopy, which requires the tracking of 1792 paths. Bertini performed this computation in 
serial in 120 seconds yielding the same set R of 36 real critical points, as required by theory |35) . 

3.3 A cubic-centered 12-bar linkage 

Consider the 12-bar spherical linkage obtained by locking the scissors of the collapsible cube 
with 12 scissors linkages presented in [52j, which is displayed in Figure 3 of [51]. Following the 
setup in [5T], we will consider the cube with side length 2 where we fix the center at the origin 
and two adjacent vertices, say P7 = (—1, 1, —1) and Ps = (—1, —1, —1). Let Pi, ■ ■ ■ ,Pe denote 
the position of the other 6 vertices yielding 18 variables. The constraints on these vertices is that 
they must maintain their initial relative distances yielding a polynomial system / consisting of 
the following 17 polynomials: 

g,, = |P, - P,f - 4, 

e {(1, 2), (3, 4), (5, 6), (1, 5), (2, 6), (3, 7), (4, 8), (1, 3), (2, 4), (5, 7), (6, 8)}; 
/i, = |P,p-3, ze {1,2,3,4,5,6}. 

The algebraic set V(/) consists of 8, 34, and 2 irreducible components of dimension 1, 2, and 
3, respectively. Table [T] presents the degrees of these components. Let V be the union of the 



.9("3, ■ 



f{u3,...,ue,U3,...,Ue) = 



. ,M6, U3, . 
■ ■ , Ue,U3, . 
U3U3 - 1 
U4U4 — 1 
U5U5 - 1 

UqUq - 1 



11 



dimension 


degree 


# components 


3 


8 


2 




4 


2 




8 


14 


2 


12 


12 


16 


1 




20 


4 




24 


1 


1 


4 


6 


6 


2 



Table 1: Irreducible decomposition of V(/) 



one-dimensional irreducible components of V(/), which has degree 36, and W be a witness set 
for V. Let Ci, . . . , Cg denote the six irreducible curves of degree 4 contained in V, and Cy and 
Cg denote the two irreducible curves of degree 6 contained in V . The components Ci, . . . , Cg 
are self-conjugate while Cy and Cg are conjugates of each other. That is, Cj U Cs contains only 
finitely many real points which must be contained in C7 H Cg . 

We used RealPoints to compute a finite set of points containing a point on each connected 
component of Vr(/) contained in V by taking 

y = [0.142,0.319,-0.286,-0.167,0.276,0.238,0.217,-0.268,-0.089, 

- 0.198, 0.287, -0.042, -0.243, 0.119, 0.309, -0.312, 0.305, 0.162]^, 

z £ R^^, 7 € C, and a G C^^ to be random of unit length. Let H be the honiotopy defined 
by (§. 

• For Step [2j we computed S using a diagonal homotopy fi^* by computing Ar\ B where 
A = V{f~ -fz) X and 

B = V(Ao(x - y) + AiV/i(a;)^ + • • • + Xn^fnixf, "qAo + • • • + a^Xn - 1). 

Since V(/ — ^z) is a curve of degree 480 and i?£ = i? fl x C^®), where £ is a random 
line in C^*, consists of 13 points, the diagonal homotopy required tracking 480 • 13 — 6240 
paths, which yielded the 1536 points in S. A witness set for V(/ — 72) was computed 
using regeneration |25j and B^ was computed using a standard 2-homogeneous homotopy. 
Overall, this computation took 5.5 minutes in parallel. 

• For Step|3] only 1440 of the 1536 paths converged and i^(E') = E\. This computation took 
18.5 minutes in serial (26 seconds in parallel) and found that the set E\ n consists of 
283 distinct points. 

• For Step |4j the homotopy membership test found that R = n i?i n consists of 24 
points, which took 80 seconds in serial. 

The set R \ Sing(/) consists of 16 points and meets Ci for i = 1, . . . , 6. This yields that 
Ci n is also one dimensional for i = 1, . . . , 6. Additionally, two points of R lie in C7 n Cg, 
one of which is presented in Figure 3 of [5T] . Each of the other six points of R, which arose 
from 30 homotopy paths in Step |4j lies in the intersection of V with some higher-dimensional 
components of V(/). 
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4 Conclusion 



Infinitesimal deformations are widely used in real algebraic geometric algorithms. By utilizing 
homotopy continuation to model the deformation, we have demonstrated that one can obtain 
an algorithm for computing a finite set of real roots of a polynomial system containing a point 
on each connected component. In particular, this algorithm computes a finite superset of the 
isolated roots over the real numbers. This is similar to basic homotopy continuation in that one 
computes a finite superset of the isolated roots over the complex numbers. The isolated complex 
roots can be identified by, for example, using the local dimension test of [5], but a similar test 
currently does not exist over the real numbers. Nonetheless, since many of the algorithms in 
numerical algebraic geometry depend only on the ability to compute a superset of the isolated 
roots, we will investigate what other computations can be performed in numerical real algebraic 
geometry building from the algorithm presented here. 
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